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t. INTRODUCTION 

The broad-scale reductions and shifts that may be expected under climate change in the availability and quality of 
stopover habitat for long-distance migrants is an area of increasing concern for conservation biologists [5]. 
Researchers generally have taken two broad approaches to the modeling of migration behaviour to understand the 
impact of these changes on migratory bird populations. These include models based on causal processes and their 
response to environmental stimulation, “mechanistic models", or models that primarily are based on observed 
animal distribution patterns and the correlation of these patterns with environmental variables, i.e. “data driven” 
models [1][6][7], Investigators have applied the latter technique to forecast changes in migration patterns with 
changes in the environment, for example, as might be expected under climate change, by forecasting how the 
underlying environmental data layers upon which the relationships are built will change over time. The learned 
geostatstical correlations are then applied to the modified data layers.. However, this is problematic. Even if the 
projections of how the underlying data layers will change are correct, it is not evident that the statistical 
relationships will remain the same, i.e. that the animal organism may not adapt its’ behaviour to the changing 
conditions. Mechanistic models that explicitly take into account the physical, biological, and behaviour responses 
of an organism as well as the underlying changes in the landscape offer an alternative to address these 
shortcomings. 

The availability of satellite remote sensing observations at multiple spatial and temporal scales [3], coupled with 
advances in climate modeling and information technologies enable the application of the mechanistic models to 



predict how continental scale bird migration patterns may change in response to environmental change [8]. In 
earlier work, we simulated the impact of effects of wetland loss and inter-annual variability on the fitness of 
migratory shorebirds in the central fly ways of North America [7]. We demonstrated the phenotypic plasticity of 
a migratory population of Pectoral sandpipers consisting of an ensemble of 10,000 individual birds in response to 
changes in stopover locations using an individual based migration model driven by remotely sensed land surface 
data, climate data and biological field data. 

With the advent of new computing capabilities enabled by recent GPU-GP computing paradigms and commodity 
hardware [4], it now is possible to simulate both larger ensemble populations and to incorporate more realistic 
mechanistic factors into migration models. Here, we take our first steps use these tools to study the impact of 
long-term drought variability on shorebird survival [9], 


2. APPROACH 

The structure of our modeling framework is shown in the attached figure. We simulate the migration trajectories 
of individual bird “particles”, including their positions, velocities, and fuel status using a spatially explicit, 
biophysical flight model. We fly birds over a user-specified Lagrangian grid and run the model at a daily time 
step to simulate spring or fall migration. The model is object-oriented, allowing the extension or incorporation of 
multiple species or mechanisms. We use evolutionary programming techniques to train the model [2], In the 
future, we hope to include interacting species to capture potential phenological changes in movement patterns 
between predator-prey species under environmental or climate change. 

We have applied the model to population sizes of 10,000 birds, a land-use grid consisting of 100,00 grid cells 
covering North America, and a daily time step. However, typical shorebird populations range in the millions and 
finer spatial resolution would provide further insight in areas of high variability under climate or land use 
changes. The evolutionary programming technique also calls for multiple simulations of bird ensemble 
movements in order to reach stability. 

Several aspects of our individual base (agent modeling) approach lend themselves to GPU-GP computing [4]; for 
example, the vector addition of bird particle and wind velocities at each time step in the spatial grid. We are 
exploring the appropriate algorithms to distribute the computational tasks between the CPU and GPU 
architectures. However, preliminary experiments suggests an increase in the number of bird “particles” simulated 
by a factor of 100 or more and increase in grid cell resolution by a factor of 10 or more will be feasible. 




{ Species Pectoral Sandpiper 
Extends Bird } 



Simulate the migration routes, timing, 
and energy budgets of individual birds 
under dynamic weather and land 
surface conditions 



{ Iterate over: Population, daily time step } 


3. CONCLUSIONS 

To illustrate our approach we present several simulations of the migration patterns of shorebirds through the 
central North American flyways under varying periods of drought variability as experienced during the past 2000 
years [9], Analysis of the drought patterns from the historical and geological record suggest that the droughts of 
the twentieth century were eclipsed several times in the past by periods of longer duration and greater spatial 
extent. These events of extreme variability likely will occur in the future and may be enhanced by anthropogenic- 
induced increased rates of climatic change. 


4. ACKNOWLEDGEMENTS 

The research described in this paper was supported under a NASA Interdisciplinary Science research proposal 
entitled, “Forecasting the effects of wetland loss and inter-annual variability on the titness of migratory bird 
species” and, in part, by a NASA Applied Information Science Technology proposal entitled “A general 
framework and system prototypes for self-adaptive Earth predictive syslems-dynamically coupling sensor with 
Earth System models.” Experiments in the application of GPULib were supported, in part, by a SB1R Phase III 
award to Tech-X from NASA entitled, ’’Exploration of GPU Computing for Advanced Bird Migration Modeling.” 
Mr. Tim Blatt is a graduate student in the Computer Science and Electrical Engineering Department at UMBC. 





5. REFERENCES 


[1] A. H. Farmer and John Weins, “Optimal migration schedules depend on the landscape and the physical environment: a dynamic 
modeling view," ./. Avian Biology’, vol, 29, pp. 405-415, 1998. 

[2] David B. Fogel, “What is evolutionary computation?,” IEEE Spectrum, vol. pp. 26-32, 2000. 

[3] Christopher S.R. Neigh, Compton J. Tucker, and John R.G. Townshcnd. “North American vegetation dynamics observed with multi- 
resolution satellite data," Remote Sensing of Environment , vol. 1 1 2, pp. 1 749- 1 772 , 2008. 

[4] John D. Owens, et al., “A survey of general-purpose computation on graphics hardware," Computer Graphics Forum, vol. 26, pp. SO- 
US, 2007. 

[5] S.K. Skagen, “Migration stopovers and the conservation of arctic-breeding Calidridine Sandpipers,” A UK, vol. 123, pp. 313-322, 2006. 

[6] James A. Smith, "Simulation of avifauna distributions using remote sensing,” IEEE International Geoscience and Remote Sensing 
Symposium, 2004 , vol. 3, pp. 1813-1814, 20-24 Sept. 2004 

[7] James A. Smith and Jill L. Deppe, "Simulating the effects of wetland loss and inter-annual variability on the fitness of migratory bird 
species," . IEEE International Geoscience and Remote Sensing Symposium, 2008, vol.. 4, pp. 8.38-841 , 7-1 1 July 2008 

[8] James A. Smith and Jill L. Deppe, "Space-based ornithology — studying bird migration and environmental change in North America," 
in Remote Sensing for Agriculture, Ecosystems, and Hydrology X, edited by Christopher M. U. Neale, Manfred Owe. Guido D'Urso, 

Proc. ofSPIE, vol. 7104, 710402 ■ CCC code: 0277-786X/08/S18 • doi: 10. l'l 17/12.801243 

[9] Connie A. Woodhouse and Jonathan T. Overpeck, “2000 years of drought variability in the Central United States,” Bulletin of the 
American Meteorological Society, vol. 79. pp. 2693-2714, 1998. 



